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Abstract: By studying some bouncing universe models dominated by a specific class 
of hydro dynamical fluids, we show that the primordial cosmological perturbations may 
propagate smoothly through a general relativistic bounce. We also find that the purely 
adiabatic modes, although almost always fruitfully investigated in all other contexts in 
cosmology, are meaningless in the bounce or null energy condition (NEC) violation cases 
since the entropy modes can never be neglected in these situations: the adiabatic modes 
exhibit a fake divergence that is compensated in the total Bardeen gravitational potential 
by inclusion of the entropy perturbations. 
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1. Introduction 

The singularity problem in cosmology Q, which is not addressed in the framework of 
inflationary paradigm [0], (note however the recent "emergent universe" proposal ||] in 
that respect) may be solved either through a Pre-Big-Bang (PBB) type phase Q or by 
means of a bounce (5|, ||. Both cases in fact demand that the universe undergoes at least one 
transition from a collapsing to an expanding epoch, although the PBB situation demands 
that this happens in the Einstein frame in which general relativity (GR) holds true. It is 
this kind of bounce that we shall be concerned with here. 

A bouncing phase may originate from quantum gravity QJ or quantum cosmology ||, 
and have been seen to occur in some string-motivated models ||. In one such model, 
based on the brane hypothesis in a five dimensional context, namely the ekpyrotic sce- 



nario [10], [11 1, the bounce in fact does not address the singularity problem since it is 



assumed to be singular. This leads to some difficulties [11] which the model still have to 
deal with (see Ref. @). 

It is worth pointing out that in some improvements of the regular bouncing models, 
not necessarily relying on general relativity at all times, perturbations may even be made 
scale-invariant jlj|, thus moving one step towards an alternative model of the Universe, 
namely one for which there is only a short period of accelerated expansion, which happens 
during the bouncing era following a contraction phase, in contrast with the relatively long 
period characteristic of usual inflationary scenarios. 
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Bouncing models without a long period of inflation afterwards, besides solving the 
singularity and the horizon problems, can be made to avoid the trans-Planckian issue [15] 
and may also provide a causal explanation, in terms of a quantum mechanical origin, for 
the large scale structures in our Universe. However, they yet depend on some theory of 
initial conditions coming from quantum cosmology and/or string theory in order to address 
the flatness and isotropization issues, which are solved naturally in usual inflationary sce- 
narios. In short, bouncing models can somehow be viewed either as potential alternatives 
to inflation, or as complementary to it, providing help in constructing a full consistent 
cosmological model for our Universe @, [|, ||, [ji]] and addressing issues which inflationary 
models do not. 

In the present work, we are interested in pointing out the relevance of entropy fluctu- 
ations during the bounce phase, in particular at the time at which the transition from an 
epoch where the null energy condition (NEC) is valid to an epoch where it ceases to be 
valid occurs (which will be called, from now on, the NEC transition), which may affect the 
observable power spectrum of adiabatic perturbations in the present expanding phase of 
the abovementioned bouncing models, including the PBB case. 

In a previous study [16|, concentrating on pure general relativity as the theory for 
describing gravity, we showed that a bouncing universe described in terms of hydrodynam- 
ical fluids was unstable with respect to adiabatic perturbations. More recently |jl3| l, we 
managed to obtain a model in which all perturbations are finite at all times, including the 
full Bardeen potential without imposing any adiabatic condition, by considering the case 
of a two fluids bouncing Universe, one of the fluids described by a scalar field. 

The case of Ref. f|] (see also Ref. @) is but one particular case of Ref. [jjl]], called 
case 4 there, with the scale factor behaving near the bounce as an even function of the 
conformal time r], at least up to the fifth order. The bounce itself, i.e. the point at which 
the scale factor derivative with respect to time vanishes, is stable in this case, in the sense 
that the perturbations are bounded while passing through it. Instabilities appear however 
as logarithmic divergences in the second derivative of the adiabatic part of the gauge 
invariant Bardeen potential, hence also in the Einstein tensor and pressure perturbations, 
at the NEC transition point. In fact, any model presenting a NEC transition faces the 
abovementioned divergence in the adiabatic Bardeen potential. 

The aim of this paper is to show, in a general framework, that the instability of adia- 
batic perturbations at the NEC transition point does not stem from a possible ill-defined 
decomposition of the pressure perturbation into adiabatic and entropy parts, but resides 
on the fact that entropy fluctuations cannot be neglected there, even for arbitrarily large 
wavelengths. Note that in the framework of the present article, entropy modes mean modes 
of mixing entropy between the various fluid constituents of the universe, as opposed to in- 
trinsic entropy modes, which we assume are absent. Hence, any attempt to calculate the 
power spectrum of perturbations on bouncing models relying only on matching conditions 
for the adiabatic Bardeen potential through the bounce is not sufficient and may lead to 
erroneous results. The entropy fluctuations, in particular at the NEC transition point, 
must be considered because they are not negligible, in fact they are extremely important 
for the stability of the model, and they may transfer power to the adiabatic perturba- 
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tions afterwards. We also settle down general conditions in which perturbations can be 
consistently defined in a bouncing model. 



2. Bouncing Background 

Let us begin by a review of the main general aspects of bouncing backgrounds. Within our 
conventions, the FLRW metric reads 

ds 2 = a 2 (ri) (dr/ 2 - 7 ydxMx j ) , (2.1) 

with the spatial three-metric 7^ given by 

70" = (l + jx 2 ) Sij, (2.2) 

and 77 being the conformal time from which one derives the cosmic time t as the solution 
of the equation adr] = dt, with a given scale factor a{rf). In Eq. ( |2.2D , the parameter /C, 
representing the spatial curvature, can be normalized to K, = 0, ±1, and x 2 = 5ijX l x^ (units 
are such that we are setting h = c = 1). 

The stress energy tensor, source of Einstein field equations, will take the form 

TMf + pJA-rf, (2.3) 

for energy density e, pressure p and 4-velocity (fluid tangent vector) u, which can be 
expressed as = (1/o)<5q, i.e. , u v = a5 u o- Einstein equations for the metric ( |2,1| ) and 
stress energy tensor ( |2.3| ) read 

H 2 + JC = £ 2 p a 2 e, (2.4) 

and 

= H 2 -H' + K = -£ 2 p a 2 (e+p), (2.5) 

where £ p = 8irG N /3 is the Planck length (G7 N being Newton constant), 7i = a' /a the con- 
formal Hubble parameter, and a prime denotes a derivative with respect to the conformal 
time 77. At the point of NEC transition we have (e + p) = 0, i.e. (3 = 0. 

We now show that in the case where the energy density can be written as a function of 
the scale factor a(r/), one can prove that this function is even around the bounce, behaving 
as 

a(v) = «o + W + er? 4 + • • • , (2.6) 

with 7] = representing the bounce. Note that such a bounce needs a NEC transition 
point 1 . Indeed, the scale factor a{rf) of a general bouncing model can be expanded around 
the bounce as follows: 

a = a + bif n + dr] 2n+1 + erf n+2 + /r/ 2n+3 + • • • , (2.7) 



1 In fact, a NEC transition point may not be present in such models if the curvature of the spatial sections 
is positive and if 2b < ao- However, in order for a realistic model to satisfy these conditions, one needs 
either to impose a tremendous amount of fine-tuning, or to have an inflationary phase between the bounce 
and the radiation dominated phase; see Ref. llfJ for details. 
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where ao > 0, and the integer n satisfies n > 1. This means we demand that a(r/) must 
be at least C 2n+3 near the bounce. In order that Eq. ( |2.7D indeed represents a bounce, the 
otherwise arbitrary parameter b must satisfy b > 0. 
The function Ti.(rj) coming from Eq. Q2.7|) reads 



H = -1 [ 2n6a r? 2n ~ 1 + (2n + l)da r/ 2n + 2(n + l)ea if n+1 + (2n + 3)/a ?? 2n+2 
% 

-2nb 2 rf n ~ x - bd(4n + l) V 4n + ■■■], (2.8) 

while /3(t?) is 

p = — [Kal - 2n(2n - l)&a r? 2n ~ 2 - 2n(2n + l)daor? 2n_1 - 2(n + l)(2n + l)ea f? 2n 

-2(2n + 3)(n + l)/a r? 2n+1 + 2n(6n - l)6 2 ?? 4n - 2 + 86dn(l + Zn)^ 11 - 1 + ■ ■ ■} .(2.9) 

In the case where the energy density can be written as a function of a, eo = eo(a), one can 
prove that a(rj) is even and that n = 1 in Eq. fl2.7| ). This can be seen as follows: from the 
energy-momentum conservation equation we obtain ade/da = — 3(e + p), implying that p 
is also a function of a, p = p(a). As a consequence, the right-hand-side of Eq. ( |2.5[ ) implies 
that j3 can also be written as a function of a, and hence it can be Taylor expanded as 

/3(a) ~ /3(a ) + Aa(ao) (6r? 2 " + ^ 2 " +1 + •••) + (^ 2 " + ^ 2n+1 + • • • )' , (2-10) 

where the ", a" indices indicate differentiation with respect to the scale factor. Equations 
( |2.9| ) and ( [2.10|) are, respectively, the left hand side (LHS) and RHS of Eq. ( |2.5|) . Comparing 
powers of the conformal time, and identifying the coefficients, one finds that the term 
containing drj 2n ~ l in the LHS has no counterpart in the RHS. Hence, d must vanish. With 
d = 0, the term containing frj 2n+1 in Eq. fl2.9|) also has no corresponding term in Eq. ( 2,1C| ), 



so that / must vanish as well. By induction, we can prove that all coefficients of terms 
containing odd powers of r] in the expansion Eq. ( |2.7| ) must be zero. Hence, a{rf) must 
be an even function of the conformal time. 2 Also, if n > 1, the term containing r] 2n ~ 2 
in Eq. (|2.9D has no counterpart in Eq. (2.1C); this imposes that n = 1, and hence a scale 



factor that behaves quadratically in i]. 

A realistic bouncing universe, i.e. one leading smoothly to a radiation dominated phase, 
cannot be modeled by means of a single barotropic fluid with constant equation of state, 
since this would imply the single fluid in question to consist in radiation, and it is well- 
known that radiation alone cannot prevent the occurrence of a singularity forming p]] , even 
for positively curved spatial sections. One could also try to describe such a bounce with 
just one fluid by allowing a varying equation of state for this fluid, assuming the equation 
of state parameter u = p/e (with e the energy density and p the pressure) to depend 
on time in such a way that lima; = 1/3 for large times. This case, somehow implicitly 
assumed in Ref. will not be treated here since it would demand complete knowledge 

2 As an alternative argument, let us note that e{rj) — e[rj{a)] — e(a) is true if and only if the inverse of 
a(rj), r){a), exists, which is possible, in a bounce framework (i.e. if a is not a monotonic function of rj), if 
and only if a{rf) is even, as a is positive. 
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of the behavior of the intrinsic entropy of the fluid, which is model-dependent. We shall 
accordingly in what follows consider the next-to-simple case of N constant equations of 
state components, later allowing for exchanges in the various modes of perturbations they 
may produce at a given scale, thereby providing the possibility to introduce mixing entropy 
modes in a model-independent way. 

We want to concentrate on the situation relevant for the rest of the paper for which the 
matter content is described by an arbitrary number of non interacting hydrodynamical per- 
fect fluids, each with a barotropic equation of state with constant ratio between the energy 
density and pressure; this is but a special case of the one where the total energy density 
can be written as a function of the scale factor, and hence with even scale factor, bouncing 
behavior of Eq (|2.6| ), and a NEC transition point. This means that the components are 
noninteracting perfect fluids with stress-energy tensors given by 

N 

T tot =E I T> ZT = +rKu" -par, (2-11) 
i=i 

relations in which u^u^ = 1 is the same timelike vector for all fluids, in agreement with 
the symmetry assumptions leading to a FLRW Universe. In the simple case with which we 
are interested here, we demand that the fluid equations of state be fixed, namely 

Pi = uji = const. (2-12) 

The total energy density and pressure that enter Einstein equations are e = J2i e i an( l 
p = ^2iPi- Energy-momentum conservation for each fluid, V^T- = implies 

e' i + 3He i (l + u i ) = => e f = Cia' 3 ^^, (2.13) 

with Cj arbitrary constants. Note at this point that the energy density in this case is indeed 
a function of a, so that, from the arguments above, a(rj) must be even with n = 1 in the 
expansion ( |2.7| ) . This situation corresponds to the so-called case (4) of Ref . |l6[] , which we 
want to investigate in greater details below. 

Specializing to two fluids and upon using Eq. (EOJ), one gets 



"'^ = £j -(l+3«*) 



ci + c 2 o 3{wi ^ 2) 



K, (2.14) 



from which it is clear that for non-positive curvature spatial section, i.e. for K < 0, it is 
only possible to have a bounce, as a point in time for which a 1 = 0, provided one of the 
fluids has negative energy; this means that for fC < 0, one of the constants c\ or c 2 ought 
to be negative. The case with positive spatial curvature can allow both constants to be 
positive, but none can vanish if one demands radiation to be present, and therefore that 
either oj\ or U2 be equal to one third. 

As an aside, let us remark that, at the background level, a scalar field cj), with dynamics 
stemming from the action 



(V^) (W) - V{#) 



-gd 4 x, (2.15) 
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is equivalent to a perfect fluid with varying equation of state since one has 

Ps = ^'f - V(0> e s = + V(4>). (2.16) 

This is no longer true at the perturbation level, unless the potential V vanishes, which is 
then equivalent to a stiff matter fluid uj = 1: when V = 0, one has indeed oc 1/a 2 — > 
oc 1/a 6 . If the kinetic terms in Eq. ( |2 . 1 6|) are negative, and including also a radiation 



fluid, one recovers the prototypical bouncing model already discussed in Ref. [13||, which 
we shall discuss later on in Sec. ^L^. In the case for which the Universe is positively curved, 
a bounce is possible even if it is dominated by a sing le scalar field jfij, [0| [0|. 



We now turn to the perturbations in this class of models, which we wish to expand on 
either a basis of a fluid-by-fluid decomposition, or into adiabatic and entropy modes. 

3. General perturbations. 

The metric for a perturbed universe with no anisotropic stress contribution can be written, 
in full generality, in the longitudinal gauge (we use the notations of Ref. p0| ) 



ds 2 = a 2 (ri) [(1 + 2$)dr? 2 - (1 - 2$) 7 ; j dxW] , (3.1) 



where $ is the gauge invariant Bardeen potential J21J, which is unique since we aim at 
describing a situation with no anisotropic pressure. The perturbed Einstein equations 
with this metric then read 

V 2 $ - m& - 3(H 2 - £)$ = 4vrGa 2 (5e (si) , (3.2) 
$" + m& + {2H' + H 2 - £)* = 47rGa 2 5p {gi) . (3.3) 

All these quantities are defined in such a way as to be gauge invariant, as emphasized by the 
superscript "(gi)" (see Ref. We shall for now on take it for granted that all quantities 

under consideration are gauge invariant, and therefore avoid the use of the notation "(gi)", 
which will be implicit. 

As a first step, in the subsection below, we show that if the matter content of a 
background bouncing model can be described by an arbitrary number of non interacting 
fluids with constant equations of state, then the full Bardeen potential and all its derivatives 
are completely regular at all times. Then, in the following subsection, we define adiabatic 



and entropy fluctuation modes. We recover the general results of Ref. [16] that there are 
always divergences in the gauge invariant adiabatic potential around the bounce, except in 
the case with bouncing behavior of Eq. (|2.6j), and around the NEC transition, but using 
a different method: instead of approximating the equations near the singular points and 
solving them afterwards, we adopt the more accurate and rigorous method of solving the 
exact equation and then expanding the corresponding solution and its derivatives near the 
singular point (note that both methods are expected to give similar results near regular 
points f22fl ). In particular, the logarithmic divergence in the second derivative of the 
adiabatic Bardeen potential at the NEC transition is reobtained for any model in which 
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this transition occurs. We then rephrase, in more general terms, the result that there 
is a class of models with NEC transition in which the full Bardeen potential and all its 
derivatives are completely regular at all times, but whose second derivative of the adiabatic 
perturbation at the NEC transition is divergent. Finally, we complete this section by 
examining the properties of the curvature perturbation 21, |2^, |24j on uniform density 
hypersurfaces associated to adiabatic perturbations. 

3.1 The case for a bounded full Bardeen potential: Fluid by Fluid decomposi- 
tion 

We now prove that if the matter content of the model is described by an arbitrary number 
iV of non interacting fluids with constant equations of state, then the Bardeen potential 
and all its derivatives are regular at all times. Using the linearity of Eqs. ( |3.2| ) and (|3.3|) , we 
can decompose the total gravitational perturbation <1> as the sum <3? = £^ $j, and construct 
the N sets of decoupled equations 

V 2 $, - - 3 (H 2 - K) $i = 4TrGa 2 5ei, (3.4) 

+ SH&i + (2H' + H 2 - /C) $i = 4irGa 2 5 Pi , (3.5) 

for each value of i = 1, • ■ ■ ,N, i.e. for each fluid. Note at this point that the decomposition 
of <I> as the sum of functions ^ is merely a mathematical tool, without any particular 
physical meaning, to prove the regularity of the the Bardeen potential. The functions <l?j 
prove convenient as they encode all the information necessary to describe the dynamical 
system as a whole. 



Substituting bpi = uJidei in Eq. (|3.5[) , and inserting it into Eq. (|3.4| ), one then obtains 
the following decoupled equations for each <£j: 

< + m{l + uji)^ - uOiV 2 ^i + [2H 1 + in 2 - K) (1 + 3^)] $i = 0, (3.6) 

Hence, 2N initial conditions are necessary to obtain the full Bardeen potential; this can 



also be seen by counting the number of degrees of freedom and constraint equations [20]. 

Models of this type, as proven above, are symmetric around the bounce. Since Ti and 
TL' are regular everywhere, and as the Ui are constants, then the coefficients appearing 
in Eqs. (|3.6| ) are completely regular at all times, in particular around the bounce and 
around a possible point of NEC transition. Then, by Fuchs property |22]| , <l?j, and all their 
derivatives are also regular at all points. Consequently, the same applies true for the full 
Bardeen potential <1> = J2i 

Note that Eqs. ([0|) and ( |3.5| ) can be cast in the form of Eqs. ( |3.6| ) if and only if the 
fluid perturbations satisfy constant equations of state, 5pi = uJi5ei, with Ui =const., i.e. if 
the fluids have vanishing intrinsic entropy perturbations; this most crucial assumption is 
not always emphasized. In the case where the equations of state contain u% which are not 
constants, Eqs. ( |3.6| ) cannot be obtained in this way, the relationship between 5pi and 5ei 
may contain divergent coefficients, as we will see below, and the proof of the regularity of 
the Bardeen potential, if possible at all, would be, at least, much more involved. Note also 
that if the parameters uji are allowed to vary, one cannot guarantee that a{rj) is even, and 
that the expansion a{rj) = clq + brf + ery 4 + • • • is valid near the bounce. 
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3.2 Unboundedness of adiabatic perturbations 

We now investigate the more usual, although alternative in the bouncing context, descrip- 
tion of the perturbations based on an expansion into adiabatic and entropy modes. Such an 
expansion works perfectly well in all other known situation encountered in cosmology |^0] . 
and it is therefore of interest to understand whether it can be rendered meaningful during 
a bouncing phase. 

Let us assume the condition, known to be valid for most fluids, that the pressure 
depends on two parameters only, namely the energy density e, and the entropy S. This 
allows the expansion 

S P = c 2 s 5e + f5S, (3.7) 



where 

9eJ s e' 3\ Hf3 / 
is the "sound velocity" (see Ref. for a discussion of this quantity, named T in that ref- 



=£=-i(i + £i (3-8) 



erence), and r = (dp/dS) e . Plugging Eq. (3J3) into (p^), and making use of the expansion 



( |3.7| ), one recovers the usual Bardeen equation 

$" + 3W(1 + eg)*' + [2W + (H 2 -K)(l + 3c 2 ) - c 2 V 2 ] 3> = |^o 2 (f8S) . (3.9) 

From now on, we shall Fourier decompose the Bardeen potential and any other relevant 
space-dependent quantities on the basis of the eigenfunctions of the Laplace-Beltrami op- 
erator as, 

$ = ^u fe (x)$ fc (77), with (V 2 + fc 2 K = 0, (3.10) 

k 

where the eigenvalues k depend on the spatial curvature K pC| ]. In what follows, we shall 
assume that such an expansion has been done for all the quantities involved, which will 
then be subsequently identified with their Fourier modes, i.e. we shall not write the index 
k, which will be implicit, and make the replacement V 2 — ► — k 2 throughout. 

The evolution equation for the adiabatic modes is given by Eq. ( p. 9] ) in which one sets 
SS = and $ — > <3? a d, namely 

<' d + 3H(1 + c 2 )^ d + [c 2 s k 2 + 2H' + {n 2 -JC){1 + 3c 2 ) ]$ ad = 0. (3.11) 



The presence of c 2 , coming from the expansion ( |3.7D , in both Eq. (3.9) and its equivalent 
for the purely adiabatic modes Eq. (|3.11| ), together with the fact, as we will see below, that 
this quantity may diverge both at the bounce and at the NEC transition point, indicate 
some possible bad behavior of $ around these points. As we will see, this is true for 
Eq. pTTD but not for Eq. (pi). 



It is well known that Eq. ( p.ll ) can be solved by means of the following usual change 
of variable |^] 

* - 3i - HU ~ H (3 12) 

* ad -^^' 6= ^m> (3 - 12) 
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which transforms the original equation into the parametric oscillator equation 



u" + ( c 2 s k 2 - — ) u = 0. 



(3.13) 



The general solution of this last equation can be constructed iteratively in the regime for 
which cjk 2 < V u = 6"/9 to yield 

^ dr 



1 = B X 



1 - k 



2 /"dr 

02 / d(J ( C S 



f + k 4 



-Q2 I dfJ ( C S 



^ ( P A( Of 



+ B * f % [i - fc2 / T d ^M) 2 /" ^ + ^ f Mc s e? f 



+ ■ 



02 J MC S 0) J p 

(3.14) 



where B\ and -B2 are constants, although usually depending on the scale k (they are in 



principle calculated through a matching with the region where c 2 k 2 > 9"/ 6 []20|), and the 
remaining terms represented by the dots are of order 0(k e ) compared with those indicated. 
We will now reobtain the divergences we got in Ref. 113], this time using expansion 



3.2.1 The bounce 



For the bounce itself, we use expansion (2/7) without the term frj 2n+3 , which is unnecessary 
here. All the following behaviors are written up to first order. 
The possible cases are: 

(1) n > 1 and K / 0. In this case, we find the following behaviors (for details, see 
Ref. ||), 



2(2n - l)(n - 1) 



3JCr] 2 



2nb 



n = — v 2n -\ 

= 1C, 



1 



z oc 



f] 



2n-2 ' 



tt oc r] 
V u oc 

7/ z 



2n~l 



(3.15) 

(3.16) 
(3.17) 
(3.18) 
(3.19) 
(3.20) 



One can see that c 2 diverges with the same power as V u . Hence, as long as k 2 <C 1, the 
expansion ( 3.14[ ) can be applied as close to the bounce as we want. Inserting Eqs. ( 3,15| ), 
(3.19) and ( |3.16j ) into ( |3.14| ) and (3.12), one can easily find that 



$ad OC T] 



2-2n 



(3.21) 



which diverges at the bounce, as stated in Ref. 16]. This divergence already appears in 
the term of order k°, and it is independent of k. 
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(2) n > 1 and /C = 0. For the special case of a flat background, the various quantities 
needed to describe perturbations are modified as 



■\bnij 1 " 



H = (3.23) 
a 

(3 = -— 2n(2n - l)/? 2 ™" 2 , (3.24) 
a 

z oc const. (3.25) 

9 oc rf, (3.26) 

Vu oc -4 • (3-27) 



Here, c 2 diverges faster than V^. Hence, expansion ( p. 14 ) cannot be applied as we approach 



the bounce. In this situation, one is obliged to approximate the equation and find its 
solution. However, in this case, a divergence already in the adiabatic perturbation itself is 
found in Ref. @, 

$ad ocr/ 3n - 2 )/ 2 e |Q| , (3.28) 



where 

a 



= k^a /[3nb(n - l)]^" 



Note there is no divergence for k = 0. 

(3) n = 1, d ^ 0, V/C. This is the case where the second derivative of a (77) is non 
vanishing and a{rf) is not even. The relevant quantities are 

C ' = ~b{2b-Ka W (3 ' 29) 

H = —rj, (3.30) 
ao 

/? = £-—, (3.31) 
1 

z oc (3.32) 

oc 77, (3.33) 

V u oc - . (3.34) 

As in the first case, c 2 diverges with the same power as V u . Again, as long as k 2 <C 1, the 
expansion (|3.14 ) can be applied as close to the bounce as we want. Inserting Eqs. ( 3JJ ), 



(|3T9|) and ( gig ) into (|1|) and (gig ) one can easily find that the adiabatic perturbation 
is finite but its first derivative diverges logarithmically, <3?^ d oc i?2^ 2 ln(r/), and its second 



derivative diverges as <£" d oc i?2^ 2 /r/, exactly as stated in Ref. [16]. There is no divergence 
for k = 0. 

(4) n = 1, V/C, and d = 0. 
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This is the case we will examine in detail in this paper. The relevant quantities are 

2 _ 8b 2 + (Kb - 12e)a 

Cs " 3b(2b-JCa ) ' 1 j 

H = —rj, (3-36) 
a 

P = K-—, (3.37) 
a 

1 

zee- (3.38) 

oc 77, (3.39) 
K = const. (3.40) 

In this case, neither c 2 nor diverge at the bounce. Hence, Eq. (3.11) is regular around 
the bounce, and so are all its solutions 3 . 

3.2.2 The NEC transition 

In what follows, we concentrate on the point where (5 = 0, so that we shift the origin of 
time : for the rest of this section, r] = when (3 = 0, and we denote by an index quantities 
evaluated at this point. 

We now assume that the scale factor around 77 = is, again, differentiable at least up 
to third order, so that the following expansion 



0(77) = a 



1 + n V + \ {ml + K) t? 2 + ia 3 7? 3 + • • ■ 



(3.41) 



holds. 4 In this relation, 03 = a"'(0)/ao, and use has been made of a"(0)/ao = 2Hq + K, 
which is a simple rewriting of j3 = 0. 

Using the expansion ( |3.41| ), we find that 

H = H + r}{H 2 + K) + \ [a 3 - H (4H 2 + 3/C)] rj 2 + 0(t? 3 ), (3.42) 

leading to 

(3 ~ [Ho {ml + 5/C) - a 3 ] 77 + O (t? 2 ) , (3.43) 
while the sound velocity takes the form 

which manifestly diverges at the NEC violating point. 



3 Note that, in all these discussions, we are not examining situations where the constants appearing in 
the relevant quantities exhibited above cancel out exactly, as this requires some fine tuning. 

4 Apart from the fine tuned case having 26 = ao, the situation for which the point of NEC transition 
coincides with the bounce [i.e. no linear term in the expansion (3.41)] is nothing but the case (2) of the 
previous subsection, which is thus already treated. 
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We will now use Eqs. ( |3,12| ) and ( 3.14 ) to evaluate the divergences in the adiabatic 
perturbation. The relevant quantities in these equations are, to leading order 



2n 2 

°~ 2 = oof \- Ho ( m ° + 5/C ) " «3] + 0( V 2 ), (3.45) 
d no 

whose behavior, combined with the divergence in c 2 , yields 
2ag[a 3 -Wo(6Wg + 5X:)J 



and finally 



H H (JC - Hp [« 3 - 3Wq (2gg + 3/C)] 



plus terms of order C(r/ 3 ). 

Inserting these expressions into Eqs. fl3,12| ) and (3.14), one obtains, in the terms with 



coefficient B\, the quantities k 2 b2C\rf' ln(ry)/2 and — fc 4 6fcf?7 2 ln(ry)/2, whose second deriva- 
tive diverges as ln(ry). This is exactly the type of divergence obtained in Ref. |L6| ] by another 
method. Here, however, the fc-dependence of the divergences is obtained more precisely. 
There is no divergence for k = 0. Note that if K = 0, the coefficient 62 vanishes, so that 
the divergence appears only at order k 4 . We will return to this point in the last section. 

The divergence in the adiabatic perturbation presented above, which is present in 
any model with NEC transition points (including the class of the precedent subsection), 
suggests that adiabatic perturbations cannot be defined in such models. We will turn to 
this point in details in Sec. [|] by concentrating on the simplest situation involving only two 
fluids. In the meantime, let us end up the setting of the general formalism for N fluids by 
looking at the curvature perturbation. 

3.3 Curvature perturbation in the adiabatic case 

Another relevant function that can be useful for calculating the primordial spectrum of 



cosmological perturbation was introduced in Ref. [23]. It is the curvature perturbation £ 
on uniform density hypersurfaces (or its generalization for non-flat background ( BST [^] ) 

Cs 3Hr7H + * = Sr i+2 *' (3 ' 48) 



where u is the total equation of state which can be evaluated by means of Eqs. (2.4) and 
fl2l)D through 



and the second equality of Eq. ( 3.48 ) stems from this relation together with the explicit 
assumption K, = 0. 

This variable is useful in particular when it comes to describing ordinary transitions 
such as the radiation to matter domination, or the reheating at the end of inflation: on 
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large (super-Hubble, often misleadingly called superhorizon fl9(1 ) scales, ( is approximately 
constant. Indeed, in this case, it is found that p4 



C ^ — -Wd, (3.50) 

e + p 

where <5p na d = 5p — c 2 5e is the nonadiabatic part of the pressure perturbation, so that ( 
is expected to be conserved for adiabatic perturbations. In this latter case, to which we 
restrict our attention in this section, it therefore suffices to evaluate it before the transition 
to obtain its value after the transition, without prior knowledge of the detailed structure of 
the transition itself. It is immediately clear however from Eq. (3.50) that this will no longer 



hold whenever a nonnegligible amount of entropy perturbation is present, or if the NEC 
violation occurs at some stage. As we have seen above, both these conditions generically 
take place in a bouncing scenario, so that care must be taken in examining this particular 
case. In a fashion similar to Eq. ( |3.12| ), one can define 



Gd = -l/V P - (3.51) 



where v gives yet another way of obtaining the gravitational potential $ through [20] 



2 ac s k 2 



with 



ifi 1 ' 2 [3 1 



Hc s V2c s 6>' 

2 



(3.53) 



where, again, care must be taken when (3 and c~ change sign, which occurs at the NEC 



, „ — , „„^_^„- ,„ g „, 

transition. Note that Eq. ( |3.51 ) is only valid provided one considers adiabatic perturba- 



tions; if entropy perturbations were present at a nonnegligible level, then Eq. ( |3.51| ) would 
have to be modified by inclusion of an extra term, proportional to the entropy perturbation 
in question (and to z), in order to define the full £. 

It is worth noting at this point that the variable v draws its importance in the theory 
of cosmological perturbations from the fact that, in the case of a single fluid (or a scalar 
field) dominating the universe, it allows to write the total (gravitational and fluid) action 
as that of a simple scalar field with varying mass, which can then be easily quantized p0|| . 
It is from this variable, by assuming a Bunch-Davies |^6[ vacuum state for its relevant 
modes, that one sets initial conditions for the cosmological perturbation. 

The effective action derivable for v permits to cast its equation of motion in the same 
form as Eq. ( p.13 ), namely 



where now the effective potential reads 

V v = — , (3.55) 
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and the RHS of Eq. ( 3.54[ ), although vanishing in the adiabatic case to which we restrict our 



attention here, contains in principle a source term proportional to the entropy perturbation. 
The perturbation v, being commonly used in the literature as a quantum scalar field, may 
be argued to possess some amount of physical significance, although it is not directly 
observable. Note also that v cannot account for all the degrees of freedom if more than 
one fluid are acting on comparable levels; in this latter case, the variable v merely encodes 
the information on the adiabatic part of the full Bardeen potential. 

The solution can be similarly expanded as in Eq. (|3.14 ) for the variable u, provided 



the replacements u — > v and 6 — > z are done. Using Eq. (|3.52| ), this gives the curvature 
perturbation directly as 



Cad OC - ~ d 

Z 



+ C 2 I dr(c. 



\ -lr I ^ i dplrjrr 



(3.56) 

where C\ and C2 are again constants depending on scale and we have dropped higher order 
terms. 

In principle, using either u or v to propagate adiabatic perturbations through a given 
period should lead to the same gravitational adiabatic potential $ a( j, and hence to the same 
primordial spectrum, i.e. the same physical predictions. In practice however, because of 
the existence of poles in the effective potentials and "sound velocity" , there are instances 



leading to discrepancies (see Refs. [19, |?J f° r a more thorough discussion). 

The results obtained in this section suggest that the very definition of adiabatic modes 
may be unattainable in such bouncing models. In fact, as explained in Sec. [2|, one needs at 
least two fluids in order to construct a realistic bouncing model. This, in turn, demands 
that the perturbations cannot be purely adiabatic: even if one considers fluids with no 
intrinsic entropy perturbation, the relative entropy contribution must be present. The 
question is if the entropy fluctuations, although necessarily present, can be considered to 
be irrelevant all along in order for the adiabatic perturbations to be sufficient to accurately 
describe the evolution of the overall perturbations, as usually assumed to be the case on 
large scales. In order to answer this question, we examine in details in the following section 
(Sec. |j) a prototypical example of a two non interacting fluids bouncing model, both fluids 
having barotropic equations of state with constant ratio. 



4. Bouncing with two fluids 

In what follows, we establish the differential equations governing the entropy perturbations 
and the Bardeen potential dynamical evolution, and show that entropy fluctuations cannot 
be neglected at the NEC transition point, and hence that purely adiabatic perturbations 
have no meaning there. To conclude this section, i.e. in |4.2| , we restrict ourselves to the 
particular case of a flat universe filled with a negative energy stiff matter and radiation, 



or, in other words, the situation equivalent to the case presented in Ref. [13]. There we 
show explicitly that, although presenting the abovementioned divergences in the adiabatic 
fluctuations, the equations governing the Bardeen potential and entropy fluctuations, when 
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taken together without any adiabatic condition, yield a perfectly regular fourth order equa- 
tion for the full Bardeen potential, which is equal to the one obtained in Ref. 13], whose 
solutions and derivatives must be regular at any time, including at the NEC transition. 
This reinforces the idea that there is no problem in defining entropy fluctuations in such 
bouncing models, but that we cannot neglect them around the NEC transition point. We 
also study, in this particular example, the behavior of curvature perturbations. 

4.1 Adiabatic and entropy modes 

We now specialize to the N = 2 case for which many calculations, in particular the evolution 
of entropy perturbations, can be done explicitly. We now have 5e = 5e\ + e)e 2 and Sp = 
5p\ + 5p2 , from which we can derive the constraint pCfl 

(a^'.i = 4nGa 2 [ei(l + ui)Su {1)i + e 2 (l + cj 2 )6u i2)i ] , (4.1) 

which will be useful later. The sound velocity now takes the simple form 

2 _ <^l(l + j^Qfl + ^2(1 + ^2^2 2 x 

which diverges at the points of NEC violation e + p = 0. With these explicit relations, one 
can calculate the entropy contribution to the pressure fluctuation. This is 



- so 1 v (1 + cji)(1 +cj 2 )eie 2 

t5S = (wi - — / ' n ' v x x 
(1 + uJi)e\ + (1 + w 2 )e2 



6e\ 5e 2 



(l + cJi)ei (l+w 2 )e 2 



(4.3) 



for which we now seek a time evolution equation. 

The perturbed fluid tangent vector, for each fluid, has components given by duo = 
—a 2 5u° = a$ and Su^i^i = —oPjijSu^ 2 % = — adiV^), thus defining the potentials v\ and 
V2- Expanding the energy-momentum conservation V fl T- lu = for both fluid to first order 
yields 

54 + m{5ei + S Pi ) + + Pi )(V 2 Vi - W) = 0, (4.4) 
for the time component, and 

[(e; +p i )d j v i ]' + 4H(ei + pi)djVi + djbpi + ( Ci +p i )d j $ = 0, (4.5) 

for a spatial component. Defining the density contrasts through 

St = — , (4.6) 

and using the background energy momentum conservation, which implies 5'j = [5e' i + 37i(l + 
Wi)5ei]/ei, one transforms Eq. (|4.4j) into 

S'^il + Ui)^' (4.7) 

Let us emphasize at this stage the well-known fact that the adiabaticity condition 

SS = 0^s = ^—(-^ M =0 (4.8) 

W1-W2 1 + UJ2 J 
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is conserved in time in the long wavelength limit k — > since (u\ — u^s' = fc 2 (t>i — U2). This 
is one of the reasons for considering adiabatic modes, the other being that the observed 



spectrum of primordial fluctuations can be reconstructed, e.g. from CMBFAST [28], with 
initial conditions deep in the radiation era satisfying s = 0, whereas isocurvature modes 
lead to significant disagreement [29] with the data [3C], although some mixture is still 
acceptable J3l[] (the situation is similar to that of having a small component of the pertur- 
bations in the form of topological defects |32| ) . 

Projecting Eq. flj.5|) along the vector kj leads to the following dynamical equation for 
the velocity potentials vf. 

v' i + H(l-3tj i )v i + Lj i —i— + $ = 0. (4.9) 



With the definition 
written as 



and Eqs. (^ 



l+Ui 

) and ( [4.9D , upon using the constraint equation 



(3/C 



k 2 )<$> 



3 

2 L - 



m [3(1 + 



(4.10) 



which is nothing but a rewriting in a convenient way of Eq. ( |3.2| ) using Eq. (O) and 
the background Einstein equation ( |2.5D , the dynamical equation for the variable s follows 
directly. This is: 

s" + H(l - 3c 2 )s' + k 2 c 2 s = %-{k 2 - 3/C)$, (4.11) 



where the definition [33] 



w 2 ei(l + uj\) + wie 2 (l + 0J2) 



z (1 + wi)ei + (1 +w 2 )e2 

has been used. It should be noted that this function, just as the "sound velocity" c 
also singular at the NEC violating point, for which (e + p) — > 0. 



(4.12) 



is 



Eq. (4.11) is the main result of this section. Note that its source term diverges at 
the NEC transition point. It means that we cannot neglect s at this point for any small 
but finite k. Hence, adiabatic perturbations cannot be defined there. Note also that the 
divergent source is of order k 2 if K, 7^ 0, and k 4 if K, = 0, exactly the orders in which the 
divergences in the adiabatic potential appear as calculated in section ||. Finally, we want 
to emphasize from Eq. (4.11) that the adiabatic case is, in most of the usual situations not 
involving NEC violation or bounces, i.e. in standard inflationary models, the only one that 
is tractable self consistently, in particular in the flat /C = situation. This is because s ~ 
solves Eq. ( [4.11 ) at leading order in k 2 . The symmetric situation, with "purely entropic 
modes" , having $ = 0, is not self-consistent because s [or 5S in Eq. (|3.9|) 1 then sources the 
gravitational potential at the same order. 

Eqs. 



( |3.9| ) and ( 4.11 ) form a closed system for the perturbation variables s and $ 
depending only on the background functions. In fact, one can express the entropy and 
sound "velocities" c z and c s simply as 



'1 + U)\ + UJ 2 + UJ1UJ2 



n 2 + K 

(3 



1. 



(4.13) 
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and 



3 V W 



whereas the left hand side of Eq. (3.9) can be given the form 



3 -ey(?6S) 



(4.14) 



(4.15) 



All these relations permit to gather the overall system into a much simplified form involving 
only the scale factor a(r/), system which is therefore particularly well suited for a more 
specific investigation. Let us accordingly now apply these relations to the particular case 
of Ref. p3| , which is, as mentioned above, equivalent to a two fluids model consisting of 
a negative energy stiff matter and radiation. We want to verify whether, without any 
adiabatic condition, Eqs ([O]) and (4.11) are consistent with the general result that the 
full Bardeen potential and all its derivatives are bounded at all times. 



4.2 A worked-out example 

We now consider the particular case of a bounce occurring for a flat (fC = 0) universe 
filled with radiation (wi = 1/3) and some negative energy stiff matter having p s = e s < 
(i.e. , u<2 = 1). Energy conservation is valid separately for both fluids, yielding e r = c r /a 4 
and e s = — c s /a 6 , with c r and c s two positive constants. Note that this stiff matter can 
also be modeled by a negative kinetic energy free massless scalar field, whose background 
dynamics and perturbations were studied in Ref. [13]. 

The background FLRW metric has a scale factor that takes the form [|R| 




where a 2 , 



(4.16) 



Note that this form satisfies the symmetry require- 



c s /c r and t/ 2 = c s /(c 2 £ 2 p ) 
ment r\ — > — rj proven in section [21. 

The relevant quantities to be calculated, namely Ti, Ti', (3, c 2 , 
example, 



H(x) = 



c 2 s (x) 



X 



770(1 +x 2 )' 
2x 2 - 7 



Ti (x) — 2/1 > > 

7/5(1 + x z y 



3(2x 2 - 



1) 



cl(x) 



2x 2 + 1 



2x 2 



1 



2x 2 



4(x 2 + l) 2 ' 



and Cg read, in this 



(4.17) 
(4.18) 
(4.19) 



where we have set for further convenience x = r]/rjQ. 

The equivalent of Eqs ( |3.6D now read [recall that for the sake of simplicity, and to avoid 
an unnecessary proliferation of indices, we note simply <I> r and the A;— modes of $ in the 
following equations, i.e. e.g. <3? r = Q r (k,r])] 



d 2 $ r 



+ 



Ax d$ r 



dx 2 1 + x 2 dx 



+ 



fc 2 2 
T + (1 + x 2 ) 2 



(4.20) 
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and 



+ 



6x d$ s 



+ 



P + 



$s = 0, 



(4.21) 



dx 2 1 + x 2 dx \ 1 + x 2 / 

where k = krjQ, with <3? = $ r + $ s , the indices r and s corresponding to radiation and 
stiff matter, respectively. As one can see, all coefficients in these equations are completely 
regular, even at the point of NEC transition, x 2 EC = 1/2. The full Bardeen potential is 
therefore regular everywhere. 

Let us now examine the evolution of the perturbations from the more usual alternative 
point of view of Eqs. (|pD and ([Hi]). From Eqs. (|415|) , (p7|) , (pLTsl) , and (glf) , one can 
turn Eqs. (|]| and into 

2x 2 -7 r o 2 



d 2 $ 



2x(4x 2 - 5) d$ 
dx 2 (1 + x 2 )(2x 2 - 1) dx 



+ 



3(2x 2 
8 



1) (l + x 2 )(2x 2 -l) 



3(l + x 2 )(2x 2 -l)' 



and 



d 2 * 
dx 2 



4x ds 2x 2 + 1 ~ 



+ 



(4.22) 



(4.23) 



2x 2 - 1 dx 2x 2 - 1 2x 2 - 1 

The coefficients of <£, d3>/dx, s and ds/dx in these equations, as well as their source 
terms, diverge as l/(2x 2 — 1) at the NEC transition points x NEC . This may imply diver- 
gences in the solutions of these equations near the NEC transition. However, taking into 
account the entropy modes (which was not possible at all in the framework of Ref. Jl6|]) 
without any adiabatic assumption allows a more precise analysis since it permits to obtain 
the fourth order equation satisfied by all the modes of i.e. not only the adiabatic ones 
(whose meaning near NEC transition is in question). The relevant equation, derived from 
Eqs. ( {Dp and ( gjg ) is 



d 4 $ lOx d 3 $ /4~ 



dx 4 1 + x 2 dx 3 



+ 



20 



1 + x 2 / dx 2 



d 2 $ 



■ + 



6xk 2 d$ 1 
1 + x 2 dx 3 



x 2 + 1 



0. (4.24) 



All the coefficients of the above equation, as could have been expected from the analysis 
in terms of Eqs. ( f4,20| ) and ( [4.21 ), are regular everywhere. Hence, using again Fuchs 
property, $ and all its derivatives must be regular at all points, i.e. there cannot be any 
divergence nowhere, no unbounded growth of the perturbation, in accordance with the 
alternative analysis based on Eqs. ( |4.2CD and ( [4.21 ). Using Eq. ( 4.17| ), one can show that 
Eq. ( 4.24 ), as expected, is the same as the one obtained in Ref. (ll| [Eq. (35) thereof], 
where the negative energy stiff matter is described in terms of the negative energy free 
massless scalar field, namely 



+ 10W 



+ 



]k 2 + 20 (H' + 27V 



$" + Wik 2 & + -k 2 [k 2 + 4(H' + 2H 2 )] $ 
3 



0. 



(4.25) 

In the same way one can obtain, for the sake of completeness, a similarly decoupled 
fourth order equation for the entropy s, which reads 



(IV) 



-2Hs"'+ 



-k 2 + 2(H' + 2H Z 



-nk 2 s'+k 2 



k 2 -2 (H' + 2H 2 ) 



s = 0. (4.26) 
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Hence, as Ti' and 7i are regular functions, the entropy fluctuations are also well behaved 
all along. 

It should be kept in mind that Eqs. (4.25) and ( [4.26 ) do in fact represent the same 
physical system, and the second does not provide any additional information not already 
contained in the first. The number of initial conditions, which is obtained from only one of 
them, is four, as one should expect in the particular case of two fluids [see the discussion 
below Eq. Q]. 

Let us now investigate the pure adiabatic (s = 0) perturbations modes. Their equation 
is given by Eq. ( |3.1l|) , which we solve by means of Eqs. ( 3.12|) and ( 3.14 ). Note that 

the approximation ( 3.14| ) 



although the term c 2 presents a (simple) pole at x 
makes perfect sense even around this point since the effective potential for u reads 

9" 8x 6 - 2x A + 20x 2 + 3 



(4.27) 



[ m (l + x 2 )(2x 2 -l)] 2 ' 

and hence presents two second order poles: the approximation actually becomes better as 
one gets closer to the NEC transition points. 

The relevant quantities are, to leading order, 



x 2 (2x 2 - 7) 



s 2a 2 (l + x 2 )(2x 2 - l) 2 
the following term being of order O (x — x NE 



1 

8^ 



), and 



||(l + x 2 ) (2x 2 -l) 



4a^ 



v / 2( 



1 



(X *^NEC ) 



13 
~6~ 



x — X 



NEC I 



[x - X„ 



(4.28) 



(4.29) 



where now the remaining term behaves as O 



of the form expected from Eqs. (|3.46|) and ( |3.45 ) for K, = 0, i.e. with 62 = 0. 



x — x 



NEC I 



. These relations are precisely 



The integrals in Eq. (|3.14|) can be evaluated explicitly, and read 



and 



which in turn yield 



1 dr 



; a oVo 



-+X + -X" 

x 3 



2«2, VO 



X 



2x 2 - 1 



+ arctan x 



(4.30) 
(4.31) 



dr 

w 



c 2 s 9 2 da 



18s 



{x 3 + 2 [2x 4 + 3(x 2 + 1)] arctan x + 2xln(x 2 + 1)} , (4.32) 



and 



c 2 2 dr 



T da 



— [10 - x 2 + 2x 4 + 2(1 - 2x 2 ) ln(x 2 + 1)] 



18(2x 2 

The divergence in the second derivative appears in the term 



(4.33) 



ad,dom 



B,k 4 



11 dr 

e 2 



c 2 s 9 2 da + 



j-2\/ 



c 2 s 9 2 dr 



T da 
I 2 



c 2 s 9 2 dp 



(4.34) 
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which, after inserting Eqs ( [4.28 ), ( 4.29 ), ( 4,32 ) and ( 4.18 ) expressed as Taylor series around 



x X* 



yields 



"NEC 



), 



(4.35) 



$ ad,dom « #lfc 4 ln (X 

as expected from the results of section [3L 

Note that in order to get Eq. ( |4.35| ), it is necessary to keep the subleading part in 
the second term of Eq. ( [4.34 ) since the leading orders of both terms are oc (x — £ NEC ) _1 
and exactly cancel each others. This result, together with the regular equations ( 4.24| ) 
and ( 4.26j ) obtained above, show that the Bardeen potential and the entropy perturbations 
given as solutions of Eqs. ( 4.22 ) and ( [4,23| ) make sense at the NEC transition, but that 
the adiabatic modes in themselves do not. However, since the divergence itself arises at 
the fourth order in k, and in the second derivative of the adiabatic part of the Bardeen 
potential, it can be argued that, in the flat case tC = at least, this adiabatic part of 
the Bardeen potential and its first derivative could be used in a consistent way to produce 
matching conditions in a bouncing scenario. 

To be complete, let us now compute the curvature perturbation as discussed in Sec. |3.3| . 
We first evaluate the background function z in the background (4.16). This gives 



^(2x 2 



1) 



3 1 



nl/2 



2x 2 - 7 



which in turns yields the following complicated form for the potential V v 



V„ 



Ux 6 + 108x 4 + 273x 2 + 98 
[rjox (x 2 + 1) {2x 2 - 7)] 2 : 



(4.36) 



(4.37) 



and thus exhibits two second-order poles whose origin can be traced back to the time 
at which the "sound velocity" c s vanishes, and another pole at the bounce itself. As is 
apparent from the solution (3.56), however, the points at which c s —* do not generate any 
divergences, and are in fact perfectly regular as far as the solution is concerned; therefore, 
we shall not consider them further. 

The solution given by Eq. ( |3.56j ) involves essentially the same integrals as those in 
Eq. ( 3.14| ), except for the last one (the triple integral), which we could not perform explic- 
itly. It turns out, however, that, contrary to the Bardeen potential, and because of the 
reverse order in which the integrands appear in the integrals of Eq. ( [3.56 ) with respect to 
Eq. ( |3.14j ), the curvature perturbation diverges at the NEC transition, and is thus not a 
proper quantity to propagate through this point. Since <3? a d, the expansion rate TC, and 
their first derivatives are all well-behaved at the NEC violating time, this fact was to be 
expected from the very definition of Cad, see Eq. ( [3.48 ). 

One point of interest to be mentioned here is the fact that, if one restricts attention 
to adiabatic perturbations, then both variables u and v exhibit divergences while passing 
through the bounce. However, those present in u stem from the NEC violation, necessary 
in this case, and can be understood as mere computational artifacts as they are exactly 
compensated in the calculation of <I> ac i, which thus appears to be the regular variable to 
consider in this instance (remember that the divergence in $ a< j appears only in its second 
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derivative, at order k 4 in the above example). In contrast, v is actually unbounded also at 
the bounce itself, as well as at the points for which the sound velocity vanishes, while Cad 
exhibits a pole around the NEC transition. 

This result merely means that the comoving hypersurfaces are not defined at the NEC 
violating points, which is hardly surprising if one recalls the definition of these hypersur- 



faces [ 2 1 1 : they are given by the requirement that the total stress-energy tensor satisfy 
if = 0. Setting Sp com the pressure perturbation with respect to these hypersurfaces, the 
comoving coordinates are then completed by the comoving time i com related to the cosmic 
time t through the relation ^4| 

at = ad V =(l- %^ dW (4.38) 



e+p 

It is clear that the divergence observed in the curvature perturbation can then be inter- 
preted as a bad choice of coordinates, the transformation to the comoving coordinates being 
singular at the NEC transition, and hence forbidden. It seems therefore that, at least in 
this example, propagating u, i.e. $ ac i, through the bounce and NEC transition makes more 
sense than propagating v, i.e. Ca,d- 

In the case one tries to match the collapsing phase to the expanding phase directly, 
without taking into account the detailed structure of the bounce, then neither $ nor ( is ac- 
tually continuous through that particular bounce; in fact, it is the combination (3 — <J>) /TC 



which passes continuously in this case [13]. This bears close resemblance with another sit- 



uation in which general relativity alone is used to perform the bounce IS, 19]; these results 
provide new example for the general framework set in Ref. [p7| . 

5. Conclusions 

In this paper, we have examined with greater details one particular case previously dis- 
cussed in Ref. [|l(]] in which adiabatic perturbations in the hydrodynamical framework 
present divergences in its second derivative around the NEC transition, see Eq. ( 4.35| ). 



even though, as inspection of Eqs. ( ]4.20| ) and ( (4.21| ) reveals, the full Bardeen potential and 
all its derivatives are perfectly regular. Eqs. (|4.20| ) and ( 4.21 ) constitute a set of equations 



which describe the Bardeen potential, hence the scalar part of the metric fluctuations, 
completely. An alternative set of equations is the one based on the definition of entropy 
fluctuations (essential for the definition of an adiabatic perturbation) given in Eq. (|3.7|), 
namely, Eqs. ( f4.23|) and ( 4.22 ). These equations present divergent coefficients essentially 



because Eq. fl3,7| ) contains Cg, which has a pole at the NEC transition, see Eq. ( |4.18| ). This 
fact suggests that the system (^23|)-(^2^) may not be appropriate to describe perturba- 
tions around the NEC transition. However, by combining these two equations, we have 
obtained a completely regular fourth order equation for the full Bardeen potential, where 
no divergences in their solutions appear at any order, even though the equations from 
which it was obtained presented coefficients with singular points at the NEC transition. 
Hence, the hydrodynamical treatment of the example based on Eqs. ( [4.23 ) and ( [4.22 ) is 



completely equivalent to the one in terms of fields worked out in Ref. |13[ , where the same 
fourth order equation was obtained. 
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Nevertheless, inspection of Eq. ( |4.23 ) shows that entropy fluctuations cannot be ne- 
glected, even for arbitrarily small but non vanishing values of the wavelength k, at the NEC 
transition time as long as its source term diverges there. Hence, adiabatic perturbations 
cannot be defined at this point and the divergence detected in Ref. [JO]] is not physically 
meaningful. That may be understood in much the same way as undamped resonances or 
shock waves: if one neglects the friction terms thanks to which the energy of the modes can 
be evacuated away, which is what we do by considering only adiabatic perturbation, one in- 
troduces artificial discontinuities (in the case of shocks) or even divergences. The adiabatic 
perturbation approximation, which does not take into account possible exchanges between 
the various fluids involved, and therefore does not count correctly the relevant degrees of 
freedom, can only make sense if one of the fluids clearly dominates over the others, or if 
there is only one fluid with complicated equation of state, as was considered in Ref. [jit]]. 

The case studied in section |4.2| is but an example for which the scale factor near the 
bounce can be expanded as a symmetric power series in the conformal time rj, namely as 
a(rj) = oo + brj 2 + erf + • • • , where oo, b and e are constants (it can be shown that the 
coefficients up to the fourth derivative of the scale factor are necessary to describe the 
passing of the perturbations through the bounce IS]). We have shown in section [2] that 
this is the only possible behavior of the scale factor in a model where the energy density can 
be written as a function of a. This situation comprehends the case of many non interacting 
fluids with constant equations of state, where we have proven that [see Eq. (|3.6| )1 the full 
Bardeen potential and all its derivatives are completely regular at all times. 

We have also shown in section || that for models with other behaviors near the bounce, 
the conclusions of Ref. |l6|] hold correct, namely, that, besides divergences at the NEC 
transition, there are divergences in the adiabatic perturbation in the bounce itself. As the 
matter content of these types of bounce cannot be modeled with many non interacting 
fluids with constant equations of state, one cannot prove in the way we did that the full 
Bardeen potential is finite in such cases, if it is. For case (1) of section ||, as the divergence 
in the adiabatic Bardeen potential comes in zero order in k, we can only conjecture that 
they are genuine divergences of the full Bardeen potential. For the other crises, duS they 
come in at least second order in k, it may be that a process analogous to the one described 
in section 4^ occurs, and the adiabatic treatment alone may not be meaningful also at 
the bounce itself. In order to see if this is the case, specific examples of this type must 
be constructed. In fact, nonsymmetric bounces may be physically more interesting than 
the symmetric ones. Indeed, some string/brane motivated models [ [lC)[] have been proposed 
in the literature, which, although mathematically inconsistent for various reasons fllTf , 
suggest that the turning point, i.e. the bounce itself, could be the time at which entropy is 
produced, in the form of radiation say, not in a smooth way. Seen from a four dimensional 
point of view, such an entropy production can be modeled by a nonsymmetric bouncing 
scenario, since the total energy density is then no longer a function of the scale factor. 
It is clear however that, according to the present analysis, the fate of linear perturbation 
theory in general relativity bouncing models cannot be guaranteed by examining adiabatic 
perturbations alone. Also, in the case linear perturbation theory really breaks down, it 
is possible that such models do have important second order contributions, which means 
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possibly strong backreaction effects. All these questions deserves further investigation. 
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